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RUNS  TEST  FOR  A  CIRCULAR  DISTRIBUTION 
AND  A  TABLE  OF  PROBABILITIES.* 
by 

Choo Ichiro  Asano  ** 


0.  Summary  • 


A  method  is  suggested  for  testing  whether  two  samples  observed  on 
a  circle  are  drawn  from  the  same  distribution.  The  proposed  test  is  a 
modification  of  the  well-known  Wald-Wolfowitz  runs  test  for  a  distribu¬ 
tion  on  a  straight  line.  The  primary  advantage  of  the  proposed  test  is 
that  it  minimizes  the  number  of  assumptions  on  the  theoretical  distribu¬ 
tion. 

*oduction. 


tlonv 

1 .  IntH< 


Circular  statistical  problems  arise  in  many  scientific  fields  such 
as  research  on  orientation  of  animals  (see  for  instance  Schmidt-Koenig 
(i960)),  time  period  analysis  for  biological  clocks,  and  rock  magnetism 
in  geology. 

As  Curray  (1956)  pointed  out,  for  very  large  samples  the  two- sample 

p 

problem  can  be  treated  by  the  use  of  the  x  -best  There  is,  however,  a 

great  need  for  a  test  that  can  be  used  when,  due  to  the  small  size  of  the 
2 

samples,  the  x  -best  is  not  applicable. 

In  1956  Watson  and  Williams  proposed  several  two-sample  tests  based 
on  the  von  Mises  distribution  (the  so-called  circular  normal  distribution) 
These  tests,  however,  are  all  parametric.  In  applications  there  is  not 
always  evidence  of  the  circular  normal  distribution. 


*  This  study  was  supported  in  part  by  the  contract  NSF-9968,  National 
Science  Foundation;  in  part  by  the  contract  Nonr  2249(05), (Nr  301-579), 
the  Office  of  Naval  Research,  Dept,  of  the  Navy,  with  The  Catholic  Univer¬ 
sity  of  America. 

**  This  paper  was  written  while  the  author  was  a  Research  Associate  in  the 
Statistical  Laboratory,  Department  of  Mathematics,  The  Catholic  University 
of  America,  during  the  academic  year  1962. 
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There  exist  also  some  non-parametric  two-sample  tests .  One  of  them, 
proposed  by  Kuiper  (i960),  is  a  modification  of  the  Kolmogorov-Smirnov 
test  for  a  circle.  Another  such  test  is  due  to  Watson  (1962)  who  applies 
a  similar  procedure  to  a  test  suggested  by  Smirnov.  At  present  the  useful¬ 
ness  of  these  tests  is  limited  since  tables  for  the  significant  values  of 
the  test  statistic  are  only  partly  available.  Finally,  a  modification  of 
the  Mann-Whitney- Cochran  test  for  a  circle  has  been  studied  by  Batschelet. 
The  result  is  still  unpublished. 

The  purpose  of  this  paper  is  to  propose  a  runs  test  for  a  circular 
distribution  in  the  two-sample  case  which  will  be  applicable  when  the 
sample  size  is  small  and  which  will  involve  r  minimum  number  of  assumptions 
on  the  theoretical  distribution.  This  at  first  appeared  to  be  only  a 
matter  of  extending  the  considerations  of  the  ordinary  theory  of  runs  on 
a  line  as  proposed  by  Stevens  (1939) >  Wald  and  Wolfowitz  (1940),  Mood  (1940) 
Swed  and  Eisenhart  (1943),  and  Wolfowitz  (1943).  But  it  soon  became  evi¬ 
dent  that  further  investigations  were  needed  The  rotatable  symmetry  for 
the  circle  changes  the  mathematical  treatment  essentially  Where  two 
separate  sets  of  observations  are  combined  in  all  the  various  possible 
ways  and  rotatable  symmetry  is  considered,  the  resulting  arrangements  were 
sometimes  identical  with  other  cases  already  obtained.  Hence,  we  must  re¬ 
duce  the  number  of  cases  by  the  number  of  arrangements  found  to  be  identi¬ 
cal  with  one  already  obtained . 

We  suppose  that  the  merit  and  the  properties  of  this  test  are  quite 
similar  to  those  of  the  ordinary  runs  test  studied  previously  for  points 
on  a  line.  For  practical  applications  the  runs  test  is  extremely  simple 
and  fast.  The  theoretical  treatment  of  the  runs  test  has  the  advantage 
that  no  discussion  is  necessary  to  Justify  the  independence  of  the  start¬ 
ing  point.  Research  on  the  power  of  the  circular  runs  test  has  not  yet 
begun .  * 
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The  numerical  table  of  the  distribution  function  of  the  test  sta¬ 
tistic  was  computed  on  an  IBM  1620. 

2 .  Probability  function  for  the  number  of  runs  on  a  circle. 

Given  two  sets  of  samples  on  a  circle,  we  wish  to  determine  the 
probability  function  for  the  number  of  runs  obtainable  by  combining  these 
two  sets  of  samples.  As  a  preliminary  step  we  consider  the  following 
partition  problem. 

Suppose  that  there  are  k  intervals  on  a  circle  and  that  in  each 
interval  there  are  n^  elements  of  the  first  sample,  where  the  n^s 
are  ordered  as  follows: 

(2.1)  nl  -  n2  2  : .  >  >  0  ,  for  1  <  k  £  N  , 

k 

where  2  n.  »  N  and  N  is  now  fixed.  After  the  consideration  of  all 
i-1  1 

possible  arrangements  of  such  k  partitioned  integers  on  a  circle,  we 
will  consider  how  to  combine  the  second  sample  among  these  k  partitioned 
integers  of  the  first  sample. 

Let 


(2.2)  S(k)  =  {n^  n2,  . ,  nk  |  ni  >  ni+i»  S  n±  =  N} 

denote  the  set  or  sets  of  ordered  Integers  determined  by  the  sample  size 
N  and  the  partition  size  k  j  e  g.,  for  N  =  4  we  have  =  {4}  , 

S( 2)  =  (3,1)  *  (2,2)  ,  S(3)  -  {2,1,1}  ,  S(4)  =  {1,1,1,1}  •  Thus  the 

set  k=l,2,...,N  contains  all  the  possible  partitions  for 

the  integer  N  . 


We  now  proceed  to  determine  the  number  of  rotatable  symmetries 
generated  from  arranging  these  k  partitioned  Integers  on  a  circle. 
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If  we  introduce  gj  such  that 


(2.3)  gj  =  (the  number  of  i's  f  n^  =  q,  i  =  l,2,...,k) 


where  n^<  Q<n^<N,  l<t<k, 

in  the  set  (S^J  k  =  • 

ment  in  S/.  \  ,  we  obtain  a  new  set 


for  j=l,2, . . . ,t  , 
we  can  characterize  each  element 

Then,  corresponding  to  each  ele- 
G(k)  Siven  "by 


(2.4) 


G(k)  =  (gl'g2' . 'gt} 


Furthermore,  let  the  common  divisors,  including  one,  of  such 
g1,g2,...,gt  be  d1,d2,...,dp  ,  where  we  put  d-j^  ■  l<d2< . <dp  and 

1  <  P  <  t  . 

If  we  omit  repeated  identical  arrangement,  then  the  number  of  cir 
cular  permutations  of  the  k  integers  is  given  by 


(2.5) 


<|)(n1,n2,  .  .  .  ,nk|N) 


where  Euler's  totient  function  Md^  ,  the  number  of  positive  integers 
less  than  dj^  and  prime  to  d^^  ,  is  given  by 


0(d±)  =  d1(l 


and  p^Pj,  •  •  •  •  are  the  different  prime  factors  of  di  ,  with  the  excep¬ 
tion  that  0(0)  =  0(l)  =1  .  This  corresponds  to  Barton  and  David  (1958, 


1962) . 

Similarly,  the  number  Y(n^,n2, . . . ,nk| d^,N)  of  symmetrical  cases 
di  # 

generated  by  (-£-) -rotation"-  is  obtained  in  the  following  way. 


The  notation 


-rotation  means  a  rotation  of 


d*  0 

( 360  x  -jl) 
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gl  S2 


S+ 


Let  denote  the  integers  -p, ....  ,-g— ,  respective¬ 

ly,  and  let  4>(m1,m2, .  .  .  ,mt  |  d^N)  be  given  by 


(2.7) 


1  uii. 

$(m^,mg, .  .  ...m^ld^N)  =  ^  2  <{)(d*)  (-p) !  (gr)>  > 

d* 


m. 


where  the  d*'s  denote  the  common  divisors  of  m1,m2,...,mt  ,  including 
one,  and  0(d*)  is  Euler's  totient  function,  with  the  exception  that 
0(0)  =  0(1)  =  1  . 

V 

Then  the  value  of  ¥(n^,n2»  •  •  •  ,nk  |  d^N)  is  given  by 

(2.8) 

Y(nx,n2, . .  .  ,nk|d1,N)  =  «t(m]L,m2, . .  .  ,mt  Id^N)-  ^2*  ^  ffn^ng,  .  .  .  ,nk|dj#N)  , 

where  2*  means  a  summation  over  all  multiples  of  d.j  among 
d1,d2,...,dp  .  Naturally,  as  a  result  of  (2-3),  we  obtain 

(2.9)  ^(n^ng, . .  •,nk|d1  =  1,N)  =  0  . 


Now,  on  the  basis  of  the  above  results,  let  us  consider  the  proba¬ 
bility  function  for  the  number  of  runs  obtained  by  arranging  simultaneous¬ 
ly  two  different  sets  of  observations  on  a  circle. 

Let  N^  and  Ng  be  the  sizes  of  the  first  and  second  sample,  res¬ 
pectively,  where  we  assume  N^<N2  without  any  loss  of  generality. 

First,  the  number  l(N^,  Ng)  of  all  possible  runs  observable  on 
a  circle  is  obtained  by  using  Euler's  totient  function  again  as  follows: 

N  +N  N  N 

(2.10)  l(N1,N2)  =  ^  ^(d)  (  2) 1  /  (■gf)  ’  (-jp) 1  ' 

where  the  summation  is  over  all  divisors,  including  one,  of  the  greatest 
common  divisor  of.  N^  and  N2  .  This  number  will  be  the  denominator  in 
the  determination  of  the  probabilities. 

Second,  in  order  to  give  the  probability  function  for  an  arbitrary 


-6- 


"2k  runs"  ,  l<k<N.  ,  let  us  determine  the  number  of  possibilities  that 
•  *  1  » 

2k  runs  will  be  observed.  Let  n£8^,  n^8^ , . . . ,n£8^  be  partitioned  inte¬ 
gers  for  N_  ,  8=1,2  .  Then  the  enumerators  of  both  non-symmetrical 

s 

arrangements  for  N1  and  N2  on  a  circle  are  given  by 


(2.11)  <D(n£s),4s) 


P 

2 

i=l 


,N) 


for  s 


1,2 


respectively,  where 


indicates  the  di  defined  pre¬ 


viously  for  s-th  sample. 

Hence,  the  number  of  arrangements  obtainable  by  combining  both 
non-symmetrical  circular  permutations  for  and  N2  is 


(2.12) 

Ii  =  k  T~T  (<t>(n£s),4s) 

1  s=l  1  d 


2  f (n^ 
i=l  1 


Furthermore,  where  the  symmetrical  circular  permutations  of  k  in¬ 
tegers  for  are  combined  with  the  non-symmetrical  circular  permutations 

for  N2  ,  the  number  of  possible  arrangements  is  given  by 


(2.13) 


Io  =  (  2 


d'^TT5  *(ni1M1>---'’41)l'3i(1)’N1>! 


a*1) 


^(^2)42),...,n(2)|N2)  -  ^  f(n^2),42),.  .  .„42)|d[2),N2)) 

Similarly,  reversing  the  roles  of  and  N2  ,  we  obtain 


(2.14)  I  =  2  . 

d[2)  di 


. . .,n^2'  |d[2;,N2) ) 
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Finally,  when  both  symmetrical  circular  permutations  for  and 
N2  are  combined,  the  munber  of  possible  arrangements  is  given  by 

{2.15)  Z.  Z  t(ni1),n<1) . n<l)d<l>,V 


,(D  „  (2) 

i  3 


i 


•I(n<a),42),.  ..,n(a)  |dS2),N2)  , 


where  B.C.D.  (fiy-fa)  indicates  the  grestest  common  divisor  of  jfjj 


and 


^  ' 


Thus  putting 


(2.16) 


!<%■  "a-  *>  *  f2)  j, 

<s[£jl  's$>  1 


we  obtain  the  probability  function  of  observing  2k  runs  by  arranging  si¬ 
multaneously  Nx  and  Np  observations  on  a  circle,  where  [sj®j )  indicates 

the  set  defined  by  (2.2)  for  s=l,2.  This  function  is  given  by 


(2.17)  Pr(k|  Nx,  N2)  -  I(N1,  N2,  k)  /  I(N1,  Ng)  for  1  <  k  <  ^  . 


3-  Distribution  function  and  the  table. 

Now  we  can  easily  obtain  the  distribution  function  for  the  number 
of  runs  on  a  circle  by  using  (2.17). 

If  2k  is  defined  to  be  the  number  of  runs,  then  the  probability 
of  an  arrangement  yielding  2r  or  fewer  runs  is 

(3-1)  P_(2k  <  2r)  =  l  P  {k|N, ,  Np)  . 

k=l  r  1  d 

The  following  table  has  been  prepared  for  use  in  testing  whether 
or  not  two  sets  of  observations  are  from  the  same  population.  The  Table 


gives  Pr(2k  <  2r}  to  5  decimal  places  for  <  N2  <  20 
range  of  from  2  to  20. 
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wlth  a 


TABLE  of  P(2r  <  2k)  ,  where  2r  is  the  number  of  runs. 


N2 

k 

Nj-2 

V3 

Nj-4 

Nx=5 

Nj-6 

N]“7 

Nx«8 

N1-9 

2 

1 

2 

•  3 

1. 

3 

1 

2 

3 

•  5 

1. 

•  25 

•  75 

1. 

4 

1 

2 

3 

4 

•33333 

1. 

.2 

.8 

1. 

.1 

.6 

1:9 

5 

1 

2 

3 

4 

5 

•33333 

1. 

. 14286 
.71429 
1. 

.07143 

.50000 

.92857 

1. 

.03846 

.34615 

.80769 

.96154 

1. 

6 

1 

2 

l 

l 

•  25 
l. 

.1 

.6 

1. 

. 04545 
.40909 
.86364 
1. 

.02381 

.26190 

.73810 

•97619 

1 

.01250 
.17500 
•  60000 
.92500 
.98750 
1. 

7 

1 

2 

3 

4 

1 

7 

.25 

1. 

.08333 

■58333 

1. 

.03333 

•33333 

.83333 

1. 

.01515 
•  19697 
.65152 
.95455 
1. 

. 00758 
.06897 
.50000 
.87879 
.99242 
l. 

. 00407 
.07724 
.38211 
.78862 
.97154 
•99593 
1. 

8 

1 

2 

3 

4 

I 

7 

8 

.2 

1. 

.06666 

•53333 

1. 

.02326 

.27907 

.76744 

1. 

.01010 

.15152 

.57576 

.92930 

1. 

.00461 

.08756 

.41014 

.82028 

.98157 

1. 

.00233 

.05128 

.29604 

.70396 

.94872 

.99767 

1. 

.00124 

.03218 

.21411 

.59406 

.89728 

.99010 

.99876 

1. 

9 

1 

2 

3 

4 

5 

6 

7 

8 
9 

.2 

1. 

.05263 

.47368 

1. 

.01818 

.23636 

.74545 

1. 

.0.0699 

.11889 

.51049 

.90210 

1. 

.00299 

.06269 

. 34328 
.76119 
•97015 
1. 

.00140 

.03497 

.23077 

.62238 

.91608 

.99441 

1. 

.00070 

.02028 

.15734 

.50000 

.84266 

•97972 

•99930 

1. 

.00037 

.01224 

.10612 

.41233 

.76067 

•95510 

.99667 

•99963 

1. 

Table  of  P{§r  <  2k}  (continued).  -9a_ 


V 

k 

V2 

V3 

Nj-4 

N  *5 

1 

N1-6 

N1-7 

N.j-8 

Nx-9 

10 

mi 

. 1666? 
1. 

.04545 

•45455 

1. 

.01370 
.  20548 

.69863 

1. 

.00498 

.09453 

.45274 

.87065 

1. 

.00198 

.04762 

.28571 

•7063$ 

.95635 

1. 

. 00086 
.02405 

.17869 

•53952 

.86426 

.98969 

1. 

. 00041 
.01354 
.11738 
.42133 
.78460 
.96344 
•99794 
1. 

.00021 
.00761 
.07672 
•  31859 
. 68141 
.92328 
.99238 
■99979 
1. 

11 

§ 

. 16667 
1. 

.03846 

.42308 

1. 

.01099 

.17582 

.67033 

1. 

.00366 

.07692 

.40659 

.84615 

1. 

.00137 

•03571 

.24176 

.65385 

.94231 

1. 

.00057 
•01753 
. 14480 
.48416 
.84050 

98303 

1. 

.00025 
. 00905 
. 08824 
•35219 
.72172 
.94344 
.99623 
1. 

.00012 
. 00488 
.05489 
. 25494 
. 60503 
.88509 
.98512 
.99940 
1. 

12 

i 

.  14286 
1. 

.03226 
. 38710 
1. 

.00862 

•15517 

.62931 

1. 

.00275 

.06314 

. 36538 
.81868 
1. 

.00096 
.02788 
. 20481 
.60288 
.92019 
1. 

.00041 
.01377 
•12353 
•  37667 
.77764 
.97327 
1. 

.00016 

.00626 

.06657 

.29339 

.65523 

.92058 

.99326 

n. 

.00007 

.00317 

.03938 

.20210 

.52754 

.84524 

.97542 

. 99866 
1. 

13 

1 

2 

3 

4 

I 

IS 

9 

10 

II 
12 
13 

. 14286 
1. 

.02857 

•  37143 
1. 

.00714 
.  .13571 
.60714 
1. 

.00210 

.05252 

.32983 

.79202 

1. 

.00070 

.02171 

•17577 

.56092 

.09756 

1. 

.00026 
.00955 
.09469 
. 37848 
.76161 
.96594 
1. 

.00010 
.00444 
.05212 
.25077 
. 60836 
.89443 
. 98978 
1. 

.00004 

.00217 

.02726 

.16379 

.47097 

.79863 

.96246 

.99756 

1. 

Table  of  P(2r  <  2k)  (continued) 


-9b- 


N2 

E 

N.j-10 

N^-ll 

Nj-12 

N1=  13 

10 

1 

2 

3 

4 

6 

7 

8 
9 

10 

.00011 

.00454 

.05123 

.24233 

.58560 

.87224 

.98119 

.99892 

.99989 

1. 

11 

1 

2 

3 

4 

1 

7 

8 
9 

10 

11 

.00006 
. 00274 
.03489 
.18492 
.50000 
.81508 
.96511 

•99726 

•99994 

1. 

. 00003 

!022ll 

.13491 
. 40997 
•74004 
.93651 
.99264 
.99966 
.99997 
i. 

12 

1 

.00003 

.00171 
.02382 
. 14006 

.41863 

.75864 
.94436 
•99427 
. 99980 
1. 

.00002 
.00094 
.01477 
.09772 
. 33001 
. 66999 
.90228 
.98523 
. 99906 
•.99998 
1. 

.00001 
.00053 
.00923 
.06795 
. 25568 
•59091 
.85373 
.97118 
•99727 
.99990 
•99999 
1. 

13 

1 

2 

3 

4 

\ 

7 

8 
9 

10 

11 

12 

13 

.00002 
.00111 
.01703 
.  10991 
; 36068 
•  69505 
..91796 
.98961 
.99956 
1. 

. 00001 
.00059 
.01011 
.07356 
.27346 
.  .59328 
.85981 
.97403 
.99783 
.99994 
1. 

. 00000 
.00032 
.00614 
.04977 
. 20682 
.50000 
•79318 
.95023 
.99386 
' .99968 
1. 

1 . 

.00000 
.00018 
.00381 
.03406 
.15657 
.41791 
.72281 
.91882 
.98688 
. 99898 
1. 

1. 

1. 

Table  of  P{2r  <  2k}  (continued) 


-9c- 


N2 

Nx=3 

Nx=4 

N1=5 

N^»6 

Nx=7 

N1=8 

Nx=9 

14 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

. 12500 
1.00000 

.25000 
. 35000 
1.00000 

.00581 

.12209 

.57558 

1. 

.00163 

.04412 
. 29902 
.76634 
1. 

.00051 

.01749 

.15123 

.52058 

.88837 

1. 

.00018 

.00722 

.07765 

.33586 

.72319 

.95558 

1. 

.00006 

.00299 

.03768 

.19696 

.51496 

.87718 

.98621 

1. 

0.00003 
.00149 
.02207 
.13483 
.41688 
.75421 
. 94762 
•99597 
1. 

15 

1 

2 

3 

4 

I 

9 

10 

II 
12 

ll 

15 

,12500 
i . 

.02174 
. 32609 
1. 

.00490 
.  10784 
•55392 
1. 

.00129 

.03737 

.27191 

.74098 

1. 

.00039 

.01392 

.13148 

•48337 

.87046 

1. 

.00013 
.00555 
. 06426 
.29910 

.68632 

.94465 

1. 

. 00005 
.00238 
.03273 
. 18450 
•51857 
.83647 
.97956 
1. 

.00002 
.00105 
. 01665 
11018 
•  36739 
•71039 
.93086 
.99385 
1. 

16 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

16 

.11111 

1. 

' 

.01961 

•31373 

1. 

.00408 

•09796 

.52653 

1. 

.00103 

.03199 

.24871 

.71827 

1. 

.00029 
. 01147 
.11445 
.45013 
.85173 
1. 

.00009 
.00424 
.05258 
.  28050 
.65761 
.93416 
1. 

.00003 

.00176 

.02573 

.15574 

.46731 

.81038 

*97359 

1. 

.00001 
. 00075 
■01274 
.09069 
.32454 
.66752 
•  01251 
•99125 
1. 

Table  of  P{2r  <  2k)  (continued) 


-9d 


N2 

B 

N1=ll 

N1=12 

Nx-13 

Nx=l4 

N1=15 

Nx=l6 

14 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

.00001 
.00073 
.01218 
.08572 
. 30612 
. 63690 
. 88879 
.98336 
•99911 
1. 

.00000 

.00037 

.00693 

.05505 

.22348 

.52665 

.81539 

•95975 

.99584 

.99985 

1. 

. 00000 
.00020 
.00404 
.03574 
. 16249 
.42873 

.73292 

.92526 

.98863 

.99921 

.99998 

1. 

.00000 

.00011 

.00242 
.02359 
.11887 
•34755 
.65245 
.88113 
•97641 
.99758 
.99989 
1 . 00000 
1. 

.00000 

.00006 

.00148 

.01567 

.08661 

.27817 

.57002 

.83131 

.95899 

.99448 

.99964 

.99999 

1. 

1. 

15 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

.00001 
.00049 
.00889 
. 06765 
.26156 
•57964 
.85665 
•97537 
.99845 
1. 

.00000 

.00024 

.00483 

.04158 

.18306 

.46602 

.76919 

.94243 

.99296 

.99969 

1. 

.00000 
.000i2 
.00272 
. 02605 
.12873 
. 36857 
.67661 
•89664 
.98221 
.99840 
•99995 
1. 

.00000 
.00006 
.00156 
.01653 
.09064 
. 28826 
.58470 
.83878 
.96229 
.99523 
.99972 

.99999 

1. 

.00000 
.00003 
.00092 
.01065 
. 06417 
. 22674 
.50000 
•77526 
.93583 
.98935 
.99908 
•99997 

1. 

1. 

. 00000 
. 00002 
.00055 
.OO696 
.04572 
. 17491 
.42407 
.70882 
.90261 
.98012 
•99774 
.99987 
1. 

1. 

1. 

16 

1 

2 

3 

4 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

.00000 
.00034 
. 00650 
.05330 
.22163 
.53034 
. 82425 
. 96604 

.99753 

1. 

.00000 

.00016 

.00342 

.03169 

.15023 

.41148 

.72249 

.92242 

.98907 

.99943 

1. 

.00000 

.00008 

.00185 

.01913 

.10204 
. 31487 
.61888 
.86321 
.97178 
•99713 
.99989 
1. 

.00000 

.00004 

.00103 

.01172 

.06947 

.23886 

.52117 

.79341 

.94465 

.99170 

.99940 

1. 

1. 

.  00000 
.00002 
.OOO58 
.00729 
.04750 
.18022 
.43300 
.71867 
.90828 
.98199 
.99807 
•99990 

1. 

1. 

.00000 
. 00001 
.00034 

. 00461 

.03277 

.13598 

.35717 

.64283 

.86402 

.96723 

.99539 

.99966 

1. 

1. 

1. 

.00000 

.00000 

.00020 

.00295 

.02277 

.10269 

.29296 

.56942 

.8l4o6 

.94726 

.99085 

•99911 

.99996 

1. 

1. 

1. 

Table  of  P{2r  <  2k}  (continued) 
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N2 

E 

Ni=3 

N-*4 

1 

Nl*5 

Nx=6 

as 

M 

4 

N1=8 

Nx-9 

17 

1 

2 

3 

4 

1 

7 

8 
9 

10 

11 

12 

i? 

S 

17 

.11111 

1. 

.0175^ 

. 29824 
1. 

.00351 

.08772 

.50877 

1. 

.00084 

.02757 

.22807 

•69591 

1. 

.00023 
.00934 
. 10048 
.41946 
.83413 
1. 

.00007 

.00340 

.04500 

.23917 

.61778 

.92067 

1. 

.00002 

.00132 

.02073 

.78185 

.96695 

1. 

.00001 
.00046 
.00843 
. 06421 
•  39051 
. 68055 
.  90844 

.98983 

1. 

18 

1 

2 

i 

I 
8 
9 

10 

II 
12 

13 

14 

15 

16 

17 

18 

. 10000 
1. 

.01563 

.28125 

1. 

.00299 

.08060 

.48657 

1. 

.00068 
.02392 
.  20984 
.67464 
1. 

.00018 
.00783 
.08861 
. 39181 

.81530 

1. 

.00005 
.00270 
.03807 
. 21487 
.58617 
.90806 
1. 

.00002 

.00101 

.01685 

•11593 

•39306 

•75357 

.95944 

1. 

.00001 

.00040 

.00771 
.06255 
. 25450 
.58722 
.87239 

.98442 

1. 

Table  of  P{2r  <  2k}  (continued) 
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N2 

B 

N.j-11 

Nj-12 

Nx=13 

Nj-14 

Nx-15 

Nj-16 

Nx-17 

Nj-18 

17 

1 

2 

? 

I 

7 

8 

9 

10 

II 
12 

14 

15 

16 
17 

.00000 
.00023 
.00484 
.04248 
. 18927 
.48284 
.79039 
•95515 
.99634 
1. 

.00000 
.00011 
.00245 
.02436 
. 12403 
.36323 
.67648 
.90023 
.98413 
.99905 
1. 

.00000 
.00005 
.00128 
.01419 
.08131 
. 26926 
.56461 

.82831 

.96016 

•99532 

.99980 

1. 

.00000 
.00002 
.00068 
. 00832 
.05300 
. 19600 
.45817 
•73906 
. 91461 
•97703 
.99888 
•99997 
1.  . 

.00000 
.00001 
.00038 
.00506 
.03548 
. 14501 
. 37450 
.66136 
.87650 
.97212 
. 99646 
•99978 

1. 

1. 

. 00000 
.00001 
.00021 
. 00309 
.02371 
. 10615 
. 30047 
.57808 
.82098 
.95053 
•  99175 
• 999p5 
.99997 

1. 

1. 

.00000 
.00000 
.00012 
. 00192 
. 01598 
.07781 
•  23974 
.50000 
.76026 
.92219 
. 98402 
. 99808 
•99988 
1. 

1. 

1. 

. 00000 
.00000 
.00007 
.00121 
. 01087 
. 05720 
. 19067 
.42902 
.69716 
.88784 
. 97278 
■99595 
. 99966 
•99999 

1. 

1. 

1. 

i 

j 

18 

1 

2 

3 

4 

I 

8 

9 

10 

II 
12 

13 

14 

15 

16 

17 

18 

.00000 

.00017 

.00365 

.03413 

.16208 

.43936 

.75618 

.94294 

•99480 

1. 

.00000 
.00008 
.00190 
.02009 
.10926 
. 34110 

.67231 

.86831 

.97674 

.99842 

1. 

.00000 
. 00003 
.00090 
.01063 
.06511 
. 23038 
.51367 
•79193 
.94651 
. 99289 
. 99964 
1. 

.00000 

.00002 

.00047 

.00609 

. 04150 
.16427 
.40980 
•  69919 
.90015 
.98053 
.99807 
•99993 
1. 

. 00000 
.00001 
.00025 
.00355 
.02665 
•11675 
. 32268 
. 60584 
.84180 
•95978 
. 99411 
•99957 
•99999 
1. 

. 00000 
. 00000 
.00014 
. 00210 
.01726 
.08296 
.25189 
•51735 
.77544 
. 93029 
.98660 
.99854 
•99992 

1. 

1. 

.00000 
. 00000 
.00007 
.00126 
.01122 
.05871 
•19437 
.44039 
.70687 

.89341 
.97481 
. 99640 
•99972 
•99999 

1. 

1. 

. 00000 
.00000 
. 00004 
. 00078 
.00746 
. 04226 
.15146 
•  36595 
. 63407 
.84856 
.95776 
.99256 
. 99924 
.99996 

1. 

1. 

1. 

. 00000 
.00000 
. 00003 
. 00048 
. 00498 
.03030 
.11712 
. 30445 

.56500 

.79950 

.93594 

.98658 

.99827 

.99987 

1. 

1. 

1. 

1. 

Table  of  P[2r  <  2k}  (continued). 


-9g- 


B 

V2 

Ni-3 

Nj-4 

*^=5 

Nj-6 

55 

n 

N1=8 

Nx-9 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

\l 

16 

17 

18 
19 

.10000 

1. 

.01429 

.27143 

1. 

.00260 

.07273 

.47013 

1. 

.OOO56 

.02089 

.19368 

.65443 

1. 

ImYrnn 

If-jLiVnrl 

Ilf 

.00004 

.00217 

.03241 

.19368 

.55652 

.89518 

1. 

.00001 
.00078 
.01384 
.09840 
•  35955 
.72517 
.95150 
1. 

.00000 
.00030 
.00608 
.05240 
. 22607 
•55025 
.85128 
.98029 
1. 

1 

2 

3 

4 

1 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

.09091 
1 . 

.01299 

.25974 

1. 

.00224 

.06726 

.45067 

1. 

.00047 

.01834 

.17921 

.63500 

1. 

. 00011 
.00524 
.06617 
. 37781 
.79217 
1. 

.00003 

.00176 

.02776 

.rr507 

.52861 

.88215 

1. 

.00001 

.00061 

.01139 

.08781 

.33216 

.69883 

.94317 

1. 

.00000 

.00022 

.00484 

.04413 

.20126 

.51552 

.82978 

.97568 

1. 

20 


Table  of  p{2r  <  2k}  (continued) 


-9h- 


N2 

B 

Nj-11 

Nj-12 

Nx=13 

Nj-14 

Nj-15 

N1=l6 

Nj-17 

19 

1 

2 

i 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
13 
16 

17 

18 
19 

.00000 
.00012 
.00278 
•02759 
.13923 
•  39974 
.72227 
.92961 
.99296 
1. 

.00000 

.00005 

.00131 

.01475 

•08533 

.28296 

.58881 

.85096 

•97111 

.99782 

1. 

.00000 

.00002 

.00064 

.00800 

.05221 

.19663 

.46484 

.75569 

•93130 

.98984 

.99942 

1. 

.00000 
.00001 
.00032 
.00445 
.03236 
.13655 
. 36228 
.65251 
.87422 
.97275 
.99694 
.99987 
1. 

.00000 

.00000 

.00017 

.00252 

.02015 

.09422 

.27763 

.55274 

.80493 

.94504 

.99089 

.99923 

99998 

1. 

.00000 
.00000 
.00009 
.00145 
.01267 
.06503 
.21088 
.46092 
.72832 
. 90685 
.97978 

•99746 

.99984 

1. 

1. 

.00000 
.00000 
.00005 
.00085 
. 00805 
.04502 
. 15944 
. 38012 
.64984 
•85963 
.96261 

.99382 

•99942 

•99997 

1. 

1. 

.00000 
.00000 
.00003 
.00050 
.00517 
.03129 
.12021 
•  31077 
.57278 

.80568 
•93907 
.98758 
. 99846 
•99989 

1. 

1. 

1. 

20 

1 

2 

2 

i 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

.00000 

.00010 

.00237 

.02495 

.13329 

.29325 

.65437 

.90591 

.98974 

1. 

.00000 

.00004 

.00097 
.01162 
.07121 
.25002 
.54802 
; 82474 
.96310 
.99693 
1. 

.00000 

.00002 

.00046 

.00613 
.04239 
.  16928 
.42305 
.71763 
.91400 
. 98601 
•99911 
1. 

.00000 
.00001 
.00022 
.00327 
.03251 
.12023 
•  32493 
.61004 
.84763 
.96378 
•99546 
.99978 
1. 

.00000 

.00000 

.00011 

.00181 

.01541 

.07611 

.23925 

.50435 

.76945 

.93146 

.99037 

.99872 

•99995 

1. 

.00000 
.00000 
.00006 
.00102 
.00947 
.05174 
.  16764 
.40312 
. 67784 
. 87930 
.97088 
•99585 
.99970 

.99999 

1. 

.00000 
.00000 
.00003 
.00057 
.00578 
.03445 
.13001 
. 32968 
.59590 
.82366 
.94790 
.99025 
.99894 
.99994 

1. 

1. 

.00000 

.00000 

.00002 

.00033 

.00361 

.02331 

.09551 

.26313 
.51456 
.76040 
.91684 
. 98084 
.99725 

•99977 
•  99999 
l. 

1. 

Table  of  P(2r  <  2k}  (continued). 


-91- 


N2 

k 

N.,-18 

V19 

Nj-  20 

19 

1 

2 

i 

I 

7 

8 

9 

10 

II 
12 

13 

14 

1% 

16 

17 

18 
19 

.00000 

.00000 

.00001 

.00031 

.00335 

.02186 

.09057 

.25254 

.50000 

•74746 

.90943 

.97814 

.99665 

•99969 

•99999 

1. 

1. 

1. 

.00000 
.00000 
.00001 
.00020 
.00236 
.01648 
.07330 
.14595 
. 39148 
.66429 
.86518 
.96258 
•99318 
.99923 
•99995 
1. 

1. 

1. 

1. 

1 

.00000 

.00000 

.00000 

2 

.00000 

.00000 

.00000 

3 

.00001 

.00001 

.00000 

4 

.00019 

.00012 

.00007 

6 

.00227 

.00146 

.00095 

6 

.01575 

.01085 

.00319 

7 

.06968 

.05157 

.03383 

8 

. 20740 

. 16499 

.12633 

9 

.43693 

.37293 

.31131 

20 

10 

.69596 

. 62707 

.56001 

11 

.87959 

.83501 

.78610 

12 

.96722 

.94843 

.92484 

13 

.99418 

.98915 

.98176 

14 

.99937 

.99854 

.99708 

15 

.99996 

.  99988 

•99971 

16 

1. 

1. 

•99998 

17 

1. 

1. 

1. 

18 

1. 

1. 

1. 

19 

1. 

1. 

20 

*. 

1. 

-10- 


4.  Numerical  example. 

To  illustrate  the  use  of  the  numerical  table  In  testing  for  random¬ 
ness  of  an  arrangement  under  the  null  hypothesis  that  two  sets  of  observa¬ 
tions  are  from  the  same  distribution,  let  us  consider  the  following  example. 
Watson  (1962)  studied  a  problem  related  to  the  migration  of  birds.  In  his 
data  the  measurements  are  given  only  to  the  nearest  5°,  as  follows: 

Control  group  (Nj  =  12):  50,  290,  300,  300,  305,  320,  330,  330, 

335,  340,  340,  355 

Experimental  group  (N2  -  14):  70,  155,  190,  195,  215,  235,  235, 

240,  255,  260,  290,  300,  300,  300. 

In  this  example,  unfortunately,  due  to  grouping  some  ties  occur 
between  values  from  the  two  different  sets  of  observations.  However, 
breaking  up  the  ties  we  can  give  upper  and  lower  bounds  for  the  number 
of  ties  in  favor  or  against  the  null  hypothesis. 

Following  Watson,  let  us  first  change  290,  300,  300  of  the  control 
group  into  285,  295,  295-  Then  the  number  of  runs  observed  is  6.  From 
the  Table  we  find  P(2r  <  6)  =  0.0040  such  that  the  two  samples  are  sig¬ 
nificantly  different  at  an  often-used  level  of  0.01. 

Second,  it  is  easy  to  see  that  for  various  ways  in  breaking  up  the 
ties  the  lower  bound  is  4  runs,  the  upper  bound  8  runs.  Thus  we  obtain 

Prob  (2r  <  4)  =  0.0002  <  Prob{2r  <  the  actual  number  of  runs) 

<  Prob (2r  <  8)  =  0.0357  • 

From  this  we  conclude  that  at  the  5$  level  the  two  groups  differ  signifi¬ 
cantly. 

I  am  indebted  to  Professor  E.  Batschelet  for  suggesting  to  work  in* 
the  field  of  circular  distribution,  and  to  Mr.  J.  F*  Gilroy,  S.J.,  for  his 
help  in  preparing  this  manuscript,  (The  Catholic  University  of  America) . 

I  wish  also  to  thank  Professor  G.  S.  Watson  (Johns  Hopkins  University)  for 
referring  me  to  Barton  and  David's  work  (1958,  1962). 
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